%ex3
%a)

%Sx = {0, 1, 2, 3, 4}
pX = zeros(1, 5);

N = 100000; %Numero de experiencias
p = 0.5; %Probabilidade de sair coroa
n = 4;   %numero de lancamentos por experiencia
m = 1;   %numero max de caras que pode sair


matriz = rand(n,N)<p; 
caras=sum(matriz);

total = 0;
p0=0;p1=0;p2=0;p3=0;p4=0;

    for c=1 : N                                                      
   
            switch (caras(1,c))
                case 0
                    p0=p0+1;
                case 1
                     p1=p1+1;
                case 2
                     p2=p2+1;
                case 3
                     p3=p3+1;
                case 4
                     p4=p4+1;
            end                 
                                                       
    end  


    pX(1,1) = p0/N;
    pX(1,2) = p1/N;
    pX(1,3) = p2/N;
    pX(1,4) = p3/N;
    pX(1,5) = p4/N;


figure(1);
hold on;

subplot(2,2,[1,2]);

x = 0:4;
stem(x,pX);

axis([-0.5 4.5 0 1.5]);
xlabel('x');
ylabel('pX');
title("Grafico Funcao de Probabilidade de X (Simulado)");